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I. INTRODUCTION 

There is a broad consensus that both on cosmological 
and on galactic scales substantial amounts of dark matter 
are needed to explain current observations. Recent data 
from type la supernovae, the large-scale-structure anal- 
ysis, and the anisotropies of the cosmic microwave back- 
ground (CMB) suggest that the dynamics of the present 
Universe is determined by a mixture of roughly 70% dark 
energy and roughly 30% (cold) dark matter (CDM). To 
understand the physical nature of these main ingredients 
of the cosmic substratum is one of the basic challenges in 
current cosmology. Numerous models have been worked 
out over the past few years, the most popular one still 
being the ACDM model which also plays the role of a 
reference setting for competing scenarios, such as the dif- 
ferent quintessence models, which are based on the dy- 
namics of a scalar field with a suitable potential [J, Q] • 
In some models interacting quintessence played a crucial 
role 0, E(. Furthermore, the cosmic medium has also 
been modelled as a gas with non-standard interactions 
given in terms of an effective one-particle force |5j,|6|,|7],|8(- 

Regarding dark matter (DM), the CDM paradigm of 
collisionless (yet undetected) supersymmetric particles 
is successful at the cosmological scale [j| HI], though 
some of its predictions at the galactic scale ( "cuspy" den- 
sity profiles and excess substructure) have been at odds 
with observations [llj, [12} . It is natural then to look for 
DM models with some type of phenomenological self- 
interaction 13] , but because of the unknown nature of 



DM in galactic halos, the type of self-interaction is also an 
open question. Hence, other models ( "w arm" DM) con- 
sider lighter more thermal particles [l4J and even non- 
thermal sources, such as coherent scalar fields have been 
contemplated [II E |H E3 . 

Since DM in galactic halos constitutes about 90 — 
95 % of the mass of galactic systems, while the latter 
have Newtonian characteristic velocities, the dynamics 
of these systems is usually examined within a Newto- 
nian framework in which, as a first approximation, vis- 
ible matter can be considered as test particles in the 
gravitational field of the DM halo. The standard models 
are either idealized configurations of collisionless matter 
obtained from Newtonian kinetic theory I19J. or mod- 
els based on n-body numerical simulations |2(J, J2lL |2^ 
that yield "universal" density or rotation velocity pro- 
files |23l[24ll25T,l26| . While the usual Newtonian approach 
may well be adequate on galactic scales, only a general 
relativistic setting can provide limits for the validity of 
the Newtonian approximation. In a general relativistic 
treatment in which DM provides the bulk of spacetime 
curvature, the motion of baryonic matter can be treated 
as test particle motion in the geometry generated by the 
DM ( "rotation velocities" become then velocities of cir- 
cular geodesic orbits [271 |28j). However, a relativistic 
approach will not only yield corrections to the Newto- 
nian case, but it also provides a framework which admits 
theoretical generalizations, e.g. the introduction of non- 
standard interactions in a systematic and self-consistent 
way. It is this property, which will be relevant for the 



present paper. Given the unknown nature of the halo 
matter, it seems to be of interest to consider various types 
of non-standard interactions within this matter and to 
explore the consequences of such an assumption. Matter 
with self-interaction may have a Newtonian type limit 
which differs from the standard Newtonian limit for non- 
interacting dark matter. This property is demonstrated 
here for two choices of a self-interaction. On this basis it 
becomes possible to discuss modifications of standard as- 
trophysical configurations where the focus of the present 
paper is on isothermal spheres. 

Gas models allow interactions to be described by ef- 
fective one-particle forces. To obtain candidates for po- 
tentially interesting self-interactions we follow a strat- 
egy similar to the one which was previously used in 
a cosmological context 0, S 0- These investigations 
demonstrated that the admissible structure of the rel- 
evant forces is severely restricted by the symmetries of 
the problem (which are the symmetries of the cosmo- 
logical principle in J5|, |y, 01) In the present paper we 
apply this formalism for the first time to a galactic scale 
in which astrophysical systems (DM halos) are necessar- 
ily inhomogcncous, hence we consider a spherically sym- 
metric and static geometry. The latter is then combined 
with the equilibrium conditions of a Maxwell-Boltzmann 
(MB) gas in the presence of a self-interacting force field. 
On the fluid level this force generates effective pressures, 
representing additional degrees of freedom that extend 
the classical Tolman-Oppenheimer-Volkov (TOV) anal- 
ysis. We demonstrate that the resulting configurations 
provide modifications of the usual isothermal sphere (the 
Newtonian limit of a "standard" MB gas), which have 
the potential to provide adequate phenomenological and 
theoretical descriptions of galactic DM. As a particular 
feature, the interaction, introduced here in a general rel- 
ativistic framework, takes the form of an additional New- 
tonian type force in the corresponding limit. This indi- 
cates possible relations to Modified Newtonian Dynamics 
(MOND) although, different from the latter, the in- 
teraction here is an interaction within the DM. 

The present paper also provides a generalization of pre- 
vious well established work on relativistic kinetic theory 
of collisionless matter applied to globular clusters J3(j 
(though our focus here is on DM halos). 

The paper is organized as follows: sections II and 
III present the general formalism of relativistic kinetic 
theory with self interactions modelled as a one particle 
force. In sections IV and V we obtain the full general 
relativistic equilibrium equations for a static, spherically 
symmetric spacetime. We specialize in section VI to a 
Lorentz-like force, linear in the particle momenta and 
without anisotropic stresses. In section VII we consider 
two possible cases of covariant couplings of this Lorentz- 
like force: firstly, to the 4-acceleration and the 4-velocity 
and secondly, to the Riemann tensor. Section VIII inves- 
tigates two types of Newtonian type limits: a "standard" 
Newtonian limit that yields modified isothermal configu- 
rations and a Newtonian type limit in which the general- 



ized force acts as an effective cosmological constant. Two 
numerical examples of configurations constructed with 
each type of Newtonian limit are discussed in section IX, 
while section X provides a summary and conclusion. 



II. KINETIC THEORY 

We assume that the particles of a collisionless relativis- 
tic gas move under the influence of a 4- force field F a . The 
equations of motion of the gas particles are J21| 



P 



dx a 
"dr ' 



F" 



Dp a 
~dV : 



(1) 



where r is a parameter along the particle worldlincs. 
Since the particle 4-momenta are normalized according to 
p a p a = —m?c 2 with constant mass m, the force F a must 
satisfy p a F a — 0. The corresponding equation for the 
invariant one-particle distribution function / = f(x,p) 
may be written as 



dx a dp a dp a 







(2) 



We shall restrict ourselves to the class of forces which 
admit solutions of J5J that are of the type of the Jiittncr 
distribution function 



/° = cxp(a + f3 aP a ), 



(3) 



where a — a(x) and (3 a is a timelike four vector. The 
particle number flow 4- vector, N a , and the equilibrium 
energy-momentum tensor, Tgq, are then defined in a 
standard way (see, e.g., [3l]]) as 



N a = c 



JdPp a f°(x,p), (4a) 



% = cJdPp a p b f(x lP ). 



(4b) 



Here, the integrals are over the mass shell m = const in 
momentum space. For the balance equations we find 

-V"„ = -mr I ,ip d({ F a ) = 0j ( 5 ) 



?eq ;b =-mc J dPp 



mc / dPf°F a 



dp a 

ad(f°F b ) 



dp b 



(6) 



where the semicolon denotes the covariant derivative. 



III. SYMMETRY CONSIDERATIONS 

The quantities V and T^q ma y be split with respect 
to the unique 4-velocity u a as 



N a = nu a , T|* ■- ^u a u b + Ph ab , (7) 



E 

7 2 



, a b 



^ab 



where h is the spatial projection tensor h a " = g" 
u a u b /c 2 . The scalars n, E, P can be, respectively, iden- 
tified with the particle number and matter-energy den- 
sities and the equilibrium pressure. Inserting yj into (J2J 
yields the constraint 



P a a ta + f3 ia . b) p a p b 



(3 a F a 



dF a 



QpO, 



(8) 



In order to make further progress, we have to introduce 
assumptions about F a . Guided by the analogy to the 
Lorentz force, we shall consider a force that is linear in 
the particle momentum, i.e. 



F a (x,p) = A ab (x)p b , 



(9) 



where A ab (x) is a function of the spacetime coordinates 
but does not depend on the particle momentum p a . Since 
the force has to satisfy p a F a = 0, it follows that, anal- 
ogously to the Lorentz-force, A ab = —A ba is valid. As a 
consequence we have dF a /dp a = and the condition © 
decomposes into 



mf3 c A c 



1 



P(a-b) = T,£{i9ab = . (10) 



As opposed to the case of the Lorentz force, the quan- 
tity A ab {x) is not specified here. We only require that 
it should be compatible with the integrability condition 
&,a-,b = &,b-,a which holds for spherically symmetric static 
spacetimes to be discussed below. The second relation 
in iJTUJ ! implies that (3 a = u a jhsT is a timelike Killing 
vector, which gives rise to the relations 



9 = < = T = , 
where T is the equilibrium temperature and 



T 



= 



(11) 



(12) 



where V a T = h a T. tb . Because of the antisymmetry of 
A ab (x) the first of the relations IJTUI) splits into 



which is completely general and quite independent of the 
action of a force. Obviously, the consistency between the 
last relation and the momentum balance (|14|t is guaran- 
teed through (|T2j) . 

The energy-momentum tensor T£k is not conserved 
and therefore it is not a suitable quantity in Einstein's 
field equations. Following earlier considerations in a cos- 
mological context (a, |6|, |7J ) , we try to map the "source" 
terms in the balances for Tgq on imperfect fluid degrees 
of freedom of a conserved energy-momentum tensor of 
the form 



rpab rpab , ab 



•-cq 



,_ab 



(17a) 

n h ab + Y, ab + 2 q^u^ , (17b) 
with u a q a -. u a Z ab = 0, £<\ = 0,(17c) 

where we can identify II as the scalar off-equilibrium pres- 
sure, E ab as the traceless anisotropic stress tensor and q a 
as the energy flux vector. Since we must have T ab -b = 0, 

then 



ab rr^ab 

t -b — -Jeq;b 



(18) 



which in the present case without an energy flux reduces 
to 



a ^bc 



h a b r 



v a n + n u a + ii b Y, ab + h" c s ca ;b 

-PV a a. (19) 



All these relations do not depend on a specific choice of 
the antisymmetric quantity A ab in J5J . Consequently, on 
the fluid level the additional interaction manifests itself 
through the appearance of the additional pressure con- 
tributions n and E ah . 

In the following we are interested in the limit in which 
the "equilibrium" state variables introduced in J7J sat- 
isfy the equation of state of a non-rclativistic Maxwcll- 
Boltzmann (MB) gas [3l| with rest mass-energy density 
P, 



a = 



V Q a = 



kuT 



' h m A 



(13) 



Using the force lO in the energy- momentum balance JHJ , 
we obtain 

(E + P) u n + V n P = nmu c K™ A cm = PV„ a . (14) 

The projection of © in direction of u a is satisfied iden- 
tically, since E = 6 = and the right hand side of the 
resulting equation also vanishes. 

It is interesting to consider the Gibbs-Duhcm relation 



dP = (E + P 



AT 



nTd 



(15) 



where [i is the chemical potential. With the identification 
igSJ this yields 



V n P = (E + P) 



T 



PV„Q , 



(16) 



with 



where 



E =p 



\/2ir zh 



> + h 



exp(a — z) 



P 



P 

z 



k B T 



(20) 



(21) 



(22) 



is the "relativistic coldness" parameter. Therefore, a is 
the fugacity scalar that is related to the chemical poten- 
tial // by 



//, 



knT 



fiz 
mc 2 



(23) 



IV. SPHERICALLY SYMMETRIC AND STATIC 
SPACETIME 

We now apply the formalism outlined so far to a spher- 
ically symmetric and static spacetime, whose metric can 
be given as 



ds 2 = - cxp (2$/c 2 ) c 2 dt 2 
+r 2 (d9 2 + sin 2 6 d<p 2 ) . 



1 - 



2GM 



dr 2 

(24) 



where <& = $(r) and M = M(r) have, respectively, units 
of velocity squared and mass. Our strategy is to com- 
bine the corresponding field equations with the equilib- 
rium conditions of the previous section. Then the 4- 
velocity and the 4-acceleration (in velocity and acceler- 
ation units) take their customary forms 



= exp(-$/ c 2 )c«5 a , 



- *'£ 



Other relevant quantities are 

cexp(-$/c 2 ) xa 



. r 



k a T 



fi a 



V 



T ~ 7 a = " Y a ' 



(25) 

(26a) 
(26b) 



where a prime denotes derivative with respect to r. Equa- 
tion (|12(l implies the well known Tolman law T ex \/—g 00 , 
which remains valid here for any choice of A ab in JjJJ . 

Particle motion in a spherically symmetric field is char- 
acterized by conservation of the of the angular momen- 
tum Ll = r 4 [{p e ) 2 + sin 2 6 (p ) 2 ]. But with F a ^ the 
particles do not follow geodesies and so y/— g 00 p° is not 
a constant of motion. In the present case the force may 
be written as F a = F° Sfi + F r <5" , hence the two nonzero 
components of F a 



A a °p b /c=(A/c)p [0 8 a r] 



F° 

F r 



-A 



Pr 



-m Ag r 



(27a) 



A P -1 



;Ae*> c2 



1 + g r 



(v r ) 



rY2 



T2 



-,1/2 



where 



A = A(r) = A r0 = -A 



(27b) 



(28) 



and we have eliminated po from g ab p a p h — —m 2 c 2 and 
used p r = m v r , as well as conservation of angular mo- 
mentum. The equilibrium condition 1)13(1 for a becomes 



3 $/ c 



:A 



- 2GM/c 2 k B T 



(29) 



Notice that A and F a have, respectively, units of accel- 
eration and force. 



V. FIELD EQUATIONS 

For the metric 1(24(1 . we have q a = and the most 
general form of the spacelike traceless tensor is given by 
E° = 0, E r r = -2E = -2S% = -2E<%, where 
E = E(r) is an arbitrary function. Thus, the momentum- 
energy tensor is 



rpti 



1 

Tr 
r 

jit/> 



<t> 



-E, 
P + U 
P + U 



2E 



(30) 



leading for ((24(1 and 1(301) to the field equations (or "equi- 
librium" TOV type equations) 



M' 



= A-K(E/c 2 )r 2 , 

M + (P + II -2E) r 3 /c 2 



G 



r(r~2GM/c 2 ) 



(31) 
(32) 



where, from l(20() and (121(1 . we know that E and P are 
functions of z and a. This has to be complemented by 
(J21 and by 



z 
z 






(33) 



$' 6E 
(n-2E)'= -a'P-(II-2£) -V-— . (34) 

c 1 r 

In the absence of the additional interaction, i.e. for II = 
E = 0, we recover the standard TOV equations while 
Eq. Q3J reduces to a ' = 0. 

Assuming A to be known, equations ((3 1(1 - 1(34(1 . comple- 
mented by ((2U|) and IJ2"T)l. constitute a system of five dif- 
ferential equations for six unknowns: M, $, z, a, II, E. 
Hence, an extra constraint is needed to render this sys- 
tem determinate. One way of achieving this determinacy 
is to provide a specific relation between II and E. The 
simplest cases are either one of 

E =0, purely isotropic case , (35a) 

II =0, purely anisotropic case (35b) 

Another possible approach follows by defining "radial" 
and "tangential" pressures as 



Pr 



11-2 2, 



P 



n + s. 



(36) 



leading to a system like l|3" H )-[|51 Jl . with J32J| and pi)l 
replaced by 



$' = G 



3 /„2 



M + (P + P r )r 6 /c 

r{r-2GM/c 2 ) '' 

/$' 2T 

P r ' = -a'P-P r ^ + — 
V c z r 



(37) 
(38) 



where 



We introduce the dimensionless variables 



r = 1-^ 



p± 
P' 



(39) 



is the "anisotropy factor" . If the latter is selected (say, 
from empirical considerations) we have also a determined 
system. 



VI. LINEAR ISOTROPIC CASE 

We consider now the case in which F a is linear in the 
momentum p a , with purely isotropic stresses character- 
ized by S = (or P r = Pj_). The system (|2TJ-(E1 
reduces to 



M' 

a' 

W 



An p 1 



2z 



= G 



M+ (pz- 1 +n/c 2 )r 3 
r(r-2GM/c 2 ) 

zre ®/ c J^ 



r-2GM/c 2 c 2 

-2 $' 



-a 



per 
z 



n 



(40a) 
(40b) 
(40c) 
(40d) 



where we used l|20[l and (|21|) to eliminate E and P in 
terms of z and p, with p = p(a, z) and z = z($) given by 



/ Z \ 3/2 

p = p c I — J exp (a - a c 

z = z c exp [($ - $ c )/c 2 ] , 



(41a) 
(41b) 



where the subscript c will denote henceforth evaluation 
at the symmetry center r — 0. It is convenient for self 
gravitating gas sources to eliminate z c in terms of the 
velocity dispersion of gas particles (19| 



k„ l r 



(42) 



and to express $ in terms of the dimensionless variable 



* 



(43) 



which reduces to the "normalized" Newtonian potential 
in the Newtonian limit. Equations l|41a|) and 141b(l then 
become 



Y EE £- 



exp 



a-a c + -C* + 7(l 



-£V\ 



,-e* 



where 



(44a) 

(44b) 

(45) 






r 



M EE 

V = 

A = 



iiTGp c _ 3 h 2 Hl p c 
M 



At: p c rl ' 

n 

Aroe^/" 



n c 

PcC 2 

2 p crlt Ae*<=/ C 



3 p c r h 2 H 2 



(46a) 
(46b) 
(46c) 
(46d) 

(46e) 



where we have used: 



p„ lt = 3h 2 H 2 /{8TrG), t: = (AnGp c /c 2 ) r 2 . 

With the definition l|42|l the parameter £ in Eq. I|45() rep- 
resents a dimensionless velocity dispersion of the parti- 
cles. A corresponding quantity, the parameter v defined 
in Eq. I)21|l . arises due to the existence of the additional 
pressure II. Both these parameters will be relevant for 
the Newtonian type limits discussed below. In terms of 
the parameters defined in l)46|) the system (|40a|l 140c 
becomes 



dM 



dx 


d* 


da; 


da 


dx 


dV 


dx 



= Y 



l + -£e 

2^ 



Oi' 



M + [£Fe e * + ^P] x 3 



x(x-2£M) 



A 



-25*, 



iv- 



x-2iM 
d* 



i re «* — 



(47a) 

(47b) 

(47c) 
(47d) 



where we have eliminated z from (|44 bj l . while Y follows 
from (|44a|l . We emphasize again, that this set of equation 
is valid for any force of the type (©■ Although for an 
explicit evaluation we might propose any functional form 
A = A(r), it is more convenient to provide a covariant 
form for A ab . 



VII. SPECIFIC FORCES 

While the assumed structure © of the force F a shares 
some features of the Lorentz force, its explicit functional 
form has been left open. The accelerations A ab are only 
restricted by the symmetry requirement A ab = — A ba . 
Since we expect this force to model an effective interac- 
tion, it may self-consistently depend on the macroscopic 
fluid quantities in addition to its (fixed, linear) depen- 
dence on the microscopic particle momentum. 



A. First case 

A simple choice for an antisymetric tensor A ab for a 
spherically symmetric fluid spacetime is 



c 



(48) 



which, for the metric (|24|) and using the variables (|46JI . 
yields 



A 



ax 



Ae«* 



M 

x 



, y ZYe** + i>V) x 



while the components of F a are 

„2 v r 
c 



F° 
F r 



-Am$'j rr e-* /c! 



f rr 



-Am$ g 



l+.9r 



(v 



r\2 



r2 



(49) 
(50a) 

-I 1/2 

(50b) 

Thus, the dimensionless quantity A is the free parame- 
ter (not necessarily constant) which determines the func- 
tional form of the force strength. 

Consequently, the set of relativistic equations to be 
solve d is (J2ZJ, with A given by (@SJ so that (|47c() and 
l|47d|l are replaced by 



da ,^d* 

ax ax 



dV 
d x 



£V + -\Y 
v 



Second case 



(51a) 

d* 

— . (51b) 

ax 



An interaction mediated by the curvature provides an- 
other choice of an antisymmetric tensor A ab : 



A ab = 



be pa 

9 ti , 



cde ' 



(52) 



where I is a characteristic length scale and R a cde are 
the components of the curvature tensor, hence the corre- 
sponding interaction is non-minimal. Now, forces which 
are proportional to the curvature tensor are usually not 
admitted in General Relativity since they violate the 
equivalence principle. Here, it is the self-consistent map- 
ping of the corresponding source terms on imperfect fluid 
degrees of freedom of a conserved energy momentum ten- 
sor according to (|17a|l - (|17c|l , which allows us to consider 
this force as an internal interaction within Einstein's the- 
ory. 

The choice (|52|) leads to 

A = e 9 tt R\ rt vru t /c 



so that the force components are F° = —Ap r /c and F r = 
Apo/c. With the help of H46[) we obtain the dimensionless 
quantity 



A = 5e«* 
(2M_ 



"+(Ytef + UP) 



Y 1 



& 



e* 



vV 



(54) 



where 



S 



4nGp C £2 = 2h^Hie p, 



— ■ (55) 



The relativistic equilibrium equations to be solved are 
then |g7J| with A given by (S3|), hence (|47c|) and (|47d|) 
must be replaced by 



ax ax 



dV 

d x 



where 



2M 



ZV + — YA 
v 



V|l + -Ce« 



d* 

da; 



vV 



(56a) 
(56b) 

(57) 



VIII. NEWTONIAN TYPE LIMITS 

Since our approach implies a non-standard degree of 
freedom, we have to define what we mean by a Newto- 
nian type limit. The free parameter £ defined in (|45(l is 
associated with the characteristic velocities of a spheri- 
cally symmetric MB gas. The second free parameter v in 
(|46|l is new and provides the velocity dispersion H/p as- 
sociated with the additional pressure II. It quantifies the 
impact of the non-standard interaction on our equilib- 
rium configuration. We expect this impact to be a small 
correction to the standard case (the case without addi- 
tional interaction), hence v <C 1. Since the dimensionless 
quantities 'J, Ai, a and V are not necessarily "small", it 
is reasonable to consider a weak field limit based on the 
conditions 



$, 



< 1 , (, < 1 and v < 1 , (58) 



9 



$' e -*/ c 



c 












V^ 


-.9oo = 


e $ c / c ^-£* 


« e*^ c 


g rr 


/$" 


<|>' 2 \ 


) + 


(9 rr Y 

2 


c 2 _ 


' V 


Qrr 


r 2t:M- 

X 


-1/2 



irrespective of how £ and v are related. The metric func- 
tions in (|24|l become up to first order in £ 



[i-e* + o(e 2 )] 

i+e^+o(a 



(53) 



(59) 



while applying (|58|l to (|44a|) we obtain 



and so the components of the force F a are 



Y = P- 



a a— Q c +* 



■o(0. 



(60) 



Newtonian type limits imply Newtonian particle veloci- 
ties v r <C c and L <C we. Under the conditions (|58|) we 
also have g rr ~ 1 and goo ~ — 1, hence the components 
of F a given by l|27[l become 



-mA, F°^F r ' 



(61) 



so that F° -C F r . Conditions (|58|) imply a weak field 
limit which can be identified with various types of "nearly 
Newtonian" conditions that arc obtained by comparing 
£ and v and by expanding all incumbent variables with 
respect to either one (or both) of these dimensionless 
ratios. 



A. The "standard" Newtonian limit 

In the general relativistic treatment of the "standard" 
MB gas (without self-interaction mediated by F a , so 
that: a — a c = v = 0) the Newtonian limit follows by ex- 
panding all relevant quantities up to first order in £. This 
leads to the "isothermal sphere" with relativistic correc- 
tions of order £. An equivalent Newtonian limit can be 
defined for the case F a ^ if together with l|58|) we have 



k > 1, 



but k ~ 0(1), (62) 



so that p c <J 2 — k II C and the relativistic correction in- 
troduced by II is of the same order but smaller than 
the thermal one, corresponding to the equilibrium hy- 
drostatic pressure P. By expanding up to O(0 we get 
modifications of the isothermal sphere characterized by 
the post-newtonian equilibrium equations 



dM 

dx 
d* 

da; 
da 

da; 
dV 

dx 
with Y given by (|60H . 



Y** +0(0, 



c 2 

M 



.v 



o(0, 



-A + O(0, 
-kY A + O(0, 



(63a) 
(63b) 
(63c) 
(63d) 



F r 



GMm 



F' 



v 



(65) 



where we have eliminated $' from l|40b|) . The Newto- 
nian structure of the force in this limit is obvious. The 
equilibrium equations are then l|6*3"| , with (j63c(l and (|63dJ) 
replaced by 



da 

dx 
dV 

da; 



A^ + O(0 



j — 

-kXY^ 

x~ 



o(0 



(66a) 
(66b) 



The fact that the additional interaction, introduced in 
a general relativistic setting, reduces to an additional 
Newtonian type force (|65|l in the appropriate limit re- 
minds of corresponding features of Modified Newtonian 
Dynamics (MOND) 29]. While the latter, however, was 
introduced as an alternative to DM, the force describes 
a self-interaction within the DM in our case. Neverthe- 
less, our result does not seem to exclude an approach 
in which the DM is replaced by baryonic matter with 
a modified Newtonian limit of a general relativistic de- 
scription, where the modification is a consequence of a 
non-geodesic motion of the (baryonic) particles of the 
medium. We plan to follow this line in future work. 



2. Second case 
For the curvature mediated force we have 

O(0, (67) 



A= SM(2M_ Y 



hence, the components of F a are 



£ 2 GMm 



2GM 



AirGp 



p0 py pr 



(68) 



The equilibrium equations are l|63|l . but now with (|63c|) 
and l|63dj) replaced by 



da 
da; 
dV_ 
dx 



S 



2M 



Y 



£ + o(0, 



(69a) 



-kSY 



2M 



Y 



£ + 0(0 ■ (69b) 



1. First case 

If we assume l|48f) , we have in the weak field limit 
M 



A = X 



o(0 



(64) 



B. Newtonian limit with a repulsive force 

If instead of l|62|) we choose (|58|l together with 

v < , \v\ = b£ > f , (70) 



then, instead of l|63|) . the equilibrium equations in the 
Newtonian limit are 



dx 


= Yx 2 +0(0, 


(71a) 


d* 

dx 


M 


(71b) 


da 
dx 


= A + 0(£), 


(71c) 


dV 
dx 


-- ZV^+{YA + 0(e), 
dx b 


(71d) 



with Y given by (JSDJ. If 1/6 < £, then equation l|71d|) 
implies V ex e«* = 1 + £* + 0(£ 2 ), hence 



n 



n c p = n c [i + o(0] 



if 



(72) 



But if 1/6 ^> £ (a more likely outcome), then the form 
of V cannot be guessed before a numerical integration. 
Still, the term \v\ V x in i)71bjl behaves as a repulsive force 
similar to a positive "cosmological constant" that can be 
associated with the matter-energy density J33] 



in r 



h r^A p crit c 



(73) 



The equilibrium equations in each case follow by inserting 
in l|71|) the appropriate form of A and taking v < 0. For 
the first case we have 



(74) 







A 


= A 


~M 

X 2 


- 


MP a; 


) 






with 






F r «-AGm 


m |n| 

r 2 c 2 




For the second case we have 


A 


= S 


x z 






1 


with 




F r f 


m 


£G 2 

z 2 


"M 


1 


11 

H 


- r 




"2M 




47Tp + 


|n|" 

c 2 . 



(75) 



(76) 



A. A modified isothermal sphere 

Let us consider a generalized force complying with 
(14811 . with the dimensionless parameter A given by 



A = 



1 



(78) 



In the "standard" Newtonian limit 1|58[) and (|62|l with 
k = 1 , we obtain hydrostatic configurations that resemble 
the isothermal sphere (case a = V = without force). 
Assuming the parameters characteristic of a galactic dark 
matter halo similar to that of the Milky Way 32]: 



<j c = 140 km/sec, 



(79) 



0.01 M Q /pc 3 



we integrate the equilibrium equations H(j3afl . I|63b|l . I|tj(ja|) 
and Ij66bl) and show the results in figure 1, where solid 
curves (marked as "MIS") denote the modified isother- 
mal sphere in comparison with the curves corresponding 
to the isothermal sphere (dotted curves marked as "IS"). 
Figure 1 (a) displays the rotation velocity profile obtained 
from (|S3, ll47hTl and ^3hl 



(77) 



V 2 r$' d* M 

% = — = -x— = J —+O{0- (80) 

crj er c dx x 

The profile V rot is qualitatively analogous to the isother- 
mal one, but reveals slight differences from the latter al- 
ready in the region occupied by visible matter (up to a 
radius of ~ 30 kpc): it is steeper near the center and 
reaches a higher maximal velocity. Instead of the "flat" 
isothermal profile, the present case shows a slight decay 
in V rot . Figure 1(b) shows the density profile, Y = p/p c , 
which reveals a wider "flat core" region with the same 
r~ 2 asymptotic decay of the isothermal sphere. Figure 
1(c) displays the function a, while figure 1(d) compares 
the curve for the pressure associated with the generalized 
force, V = n/II c with that for the hydrostatic pressure 
P/P c — Y e^*. While P ex Y, as in the isothermal case, 
V decays exponentially. 

It is worthwhile mentioning that similar models of 
modified isothermal spheres can be obtained with the 
curvature coupled force l|52|) , though we will use this type 
of force to model an effective cosmological constant. 



IX. NUMERICAL EXAMPLES 

We examine separately two configurations, one for 
each form of the generalized forces (|48|l and (|52|l . These 
examples aim at illustrating how the formalism we 
have introduced works in practice. Although these 
examples are not meant to describe "realistic" models of 
equilibrium systems, they show how generalized forces 
in collisionless gases might influence known velocity and 
density profiles and hence are of interest in astrophysical 
applications. 



B. Cosmological constant from curvature coupling 

Assuming a curvature mediated force l|52J). we exam- 
ine now the case of a weak field limit that resembles 
an isothermal sphere in equilibrium with a cosmolog- 
ical constant (i.e. a cosmological field with matter- 
energy density E A = -P A = Q x h 2 p CIit c 2 |33|). This 
suggests a generalized force of cosmic origin, character- 
ized by very large length scales. In fact, from l|55|) with 
H = 100km/(secMpc), 



S 



7.4 x 10" 



/r^ 2 



Pc 



(Mpc) 2 p otH 



0.0015 



(Mpc) 



r, (81) 
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FIG. 1: Panels (a) and (b) depict the rotation velocity and 
density profiles for the modified isothermal sphere (MIS), in 
comparison with those of the isothermal sphere (IS), for a 
Milky Way size galactic halo under the assumption of the 
coupling 148 1 with A = 1/(1 + x 2 ). Panel (c) displays the 
function a(x) and panel (d) provides a comparison between 
the hydrostatic pressure P/P c — Y e^ 9 and the pressure V = 
II/n c associated with the generalized force. See section IX. 
A. for more details. 



FIG. 2: Panels (a), (b) and (c) depict the curves (marked 
as "AIS") of the normalized potential 9, squared rotation 
velocity and density profile for the case in which the curvature 
coupled generalized force 1521 acts as an effective cosmological 
constant. The corresponding curves of the isothermal sphere 
(IS) are shown for comparison. Panel (d) displays the function 
a(x). Notice that V 2 at = at the same radius (r ~ 3.6 Mpc) 
where d^/dx = 0. See section IX. B. for more details. 
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where we have used h = 0.7 and (|79|l . it follows that 
5 = 1 corresponds to a characteristic length scale of 
25.87 Mpc, roughly an order of magnitude larger than 
the size of the largest galactic clusters in virialized equi- 
librium (it is roughly the size of a supercluster |3J|). Fol- 
lowing l|73|l . the ratio of central hydrostatic pressure to 
the pressure of the A-field is 



b = 



h 2 ^A Pent c 2 
PcC 2 c 



40 



(82) 



where we have used (J75J, as well as h — 0.7 and £!a = 0.7. 
Hence, in order to have a curvature mediated force acting 
as a cosmological constant, it is reasonable to assume 
S = 1, together with conditions (|70|l with 6 = 40. Notice 
that the rotation velocity profile is no longer given by 
JHOJ, but by 



V 2 



M 



W\Vx 2 + O(0- 



(83) 



We integrate the equilibrium equations (|71|l for the 
values ((79|) . with A given by (J7T)|) . displaying the results 
in Fig. 2. Figures 2(a)-2(c) show how the normalized 
potential, ^, the velocity squared and density profiles, 
V r 2 t /a 2 and Y, are practically indistinguishable from 
the curves of these same variables in the case of an 
isothermal sphere coupled to a cosmological constant 
(see [33j ] for a comparison and a detailed discussion 
of this system). As figure 2(a) shows, the function "f 
oscillates, indicating "turning points" where <&' = that 
separate regions with different signs of $'. The density 
profile in figure 2(c) does not decay asymptotically as 
1/r 2 , as the isothermal profile, but oscillates around 
the density value of the "background" A-field. At the 
first turning point of $' the squared velocity V 2 t (figure 
2(b)) vanishes and becomes negative (hence circular 
geodesic orbits cannot exist for larger r). Notice that 
this first turning point occurs at r « 3.6 Mpc, about an 
order of magnitude larger than the physical radius of an 
isothermal halo similar to the Milky Way. As discussed 
in |23, this fact allows one to ignore the effects of a 
cosmological constant in the equilibrium of galactic DM 
halos. A phenomenological analogue of a cosmological 
constant can also be achieved by applying conditions 
(j7U|l to a force given by the coupling P%|) . though we 
have preferred to illustrate this analogue with the case 
(|52|l because it offers a more natural association with a 
cosmic force acting on very large length scales. 



X. DISCUSSION AND CONCLUSION 

We have derived the hydrostatic equilibrium conditions 
for a gas with a specific class of self-interactions in a 
spherically symmetric spacetime. This self-interaction is 
described by a Lorentz-like force © that is linear in 



the particle momentum and characterized by a Maxwell- 
like antisymmetric tensor A ab . Within this framework, 
we have considered two possible covariant forms for this 
force that are adequate for the static and spherically sym- 
metric spacetime geometry: the one given in equation 
(|48|l where this tensor is proportional to v\ a u h \ the other 
given in equation (|52|) . coupling it to the Riemann ten- 
sor. The equilibrium equations in each case were pre- 
sented in terms of dimensionless variables, which allowed 
us to obtain the Newtonian type limit in a natural way 
by expanding with respect to the dimensionless param- 
eters £ and v, which characterize the dispersion veloc- 
ities associated with the hydrostatic ideal gas pressure, 
P = nksT, and the pressure II, respectively, the latter 
being the result of the self-interaction on the hydrody- 
namic level. We have provided two numerical examples 
of Newtonian type weak field configurations, a modified 
isothermal sphere and an example of how the curvature 
coupled force can act phenomenologically as an effective 
cosmological constant. 

While we do not claim the examples introduced in 
section IX and depicted by figures 1 and 2 to describe 
"realistic" configurations, we believe that these simplest 
possible cases nevertheless do illustrate how configura- 
tions of astrophysical interest might arise from the pre- 
sented formalism. In fact, these models, especially that 
of the modified isothermal sphere shown in figure 1, do 
convey interesting information worth discussing. For ex- 
ample, numerical simulations pol I2I I2I I2H I2I I25I l2r| 
yield halo configurations characterized by empiric den- 
sity and velocity profiles that differ from the isothermal 
profiles, at least in the inner halo region containing a disk 
of baryonic matter (which can be considered as "tracers" 
of the halo gravitational field) . The non-isothermal ve- 
locity profile associated with the Navarro-Frenk- White 
(NFW) and analogous simulations I20l l2ll |22J |. as well 
as other velocity profiles J23, Hj, [23, |2|j, are qualita- 
tively similar to the velocity profile shown by the solid 
curve in Fig. 1(a), corresponding to an MB gas with a 
generalized force of the type l|48|l . Surprisingly, the em- 
piric density profiles of these simulations show a "density 
peak" towards the halo center which does not match the 
density profile of Fig. 1(c) that shows a "flat density 
core" . Numerical simulations also yield an asymptotic 
decay p oc 1/r 3 , that is different from the isothermal de- 
cay p oc 1/r 2 of figure 1(c) (though it is difficult to verify 
the asymptotic DM density decay by means of observa- 
tions). However, observations in low surface brightness 
(LSB) galaxies |35l l3q (supposedly dominated by DM) 
do not reveal the "density peak" predicted by numerical 
simulations, but an isothermal "flat core" (this is still a 
hotly controversial issue). Hence, the simple toy model 
we have proposed in section IX shows a rotation velocity 
profile compatible with numerical simulations, but with- 
out the controversial "density peak". It is then quite 
possible that more elaborated models of a MB gas with 
self-interaction could provide a reasonable phenomeno- 
logical fit to observations, especially if we consider the 
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case with anisotropic stresses since more realistic galactic 
halos are very unlikely to correspond to isotropic config- 
urations. Further research in this direction is currently 
being undertaken. 

It is expedient to point out that the interactions stud- 
ied in this paper seem to share features with MOND [2Jj • 
This could provide the theoretical basis for an approach 
in which the DM is replaced by baryonic matter with 
a modified Newtonian limit of a general relativistic de- 
scription. In such a setting modifications of Newton's 
law would be the result of a non-geodesic motion of the 
(baryonic) particles of the medium. 

Finally, regarding the phcnomcnological modelling of a 
cosmological constant by means of the curvature coupled 
force l)52[l . we have shown that this self-interaction can 
reproduce similar effects as those of previous studies of 
a cosmological constant (a cosmic A field) in hydrostatic 
equilibrium with a MB gas (see (33). On the other hand, 



the gas dynamics with self-interactions of the type that 
we have presented in this paper has been successfully 
applied to dark energy models at the cosmological 
scale [a, Ll Ll @- Thus, the study of the interplay and 
correspondence between galactic scale models, as those 
derived here, and those describing the cosmic dynamics 
dominated by dark energy provides an interesting and 
relevant avenue for further research. 
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